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Abstract 

r» , 

i This paper deals with the approximation of the spectrum of linear 

. and nonautonomous delay differential equations through the reduction of 

' the relevant evolution semigroup from infinite to finite dimension. The 

focus is placed on classic collocation, even though the requirements that a 
^ . numerical scheme has to fulfill in order to allow for a correct approxima- 

' tion of the spectral elements are recalled. This choice, motivated by the 

analyticity of the underlying eigenfunctions, allows for a convergence of 
infinite order, as rigorously demonstrated through a priori error bounds 
when Chebyshev nodes are adopted. Fundamental applications such as 
, determination of asymptotic stability of equilibria (autonomous case) and 

5^ ■ limit cycles (periodic case) follow at once. 
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1 Introduction 

In the recent decades, Delay Differential Equations (DDEs), and more general 
functional differential equations, have attracted the attention of diverse scientific 
communities, beyond that of mathematicians, ranging from automatic control to 
physics, through population dynamics and bio-mathematics, to name a few. A 



1 



central question from a dynamical point of view is that of being able to determine 
or foresee the asymptotic stability of equilibria and periodic orbits of nonlinear 
systems through suitable linearization. Despite the great effort, the utilization 
of well-established analytical results such as spectral bounds and Lyapunov 
stability turns out to be rather lacking from the practical point of view of 
applications and, at best, suitable for restricted sub-classes (e.g. single discrete 
delay, second order systems, etc.). A major difficulty is clearly manifested in the 
infinite dimensional nature of these models arising when their time evolution 
is described in a suitable state space. As a natural consequence [26, p. 109], 
a number of approximation techniques have been proposed, mostly based on 
computing the characteristic values (roots, multipliers, Lyapunov exponents) 
associated to the system, see e.g. [8, 12, 19, 20, 22, 27, 29, 37, 39]. 
When investigating on the stability (but not only) of 



where T : [0, oo) x C — > C d is linear with C the state space of elements x t € C 
defined as 



according to the standard Hale-Krasovskii notation [26, 31], r being the max- 
imum delay and d the number of equations, the state space description of the 
model is advantageous, and the classic literature resorts to the Banach space 
of continuous functions C := C(— t, 0;C d ), [2, 18, 26, 40]. This choice seems 
to be motivated by the fact that, for rather general selections of the space of 
initial data, the "smoothing effect" [2] makes the (forward) solution be continu- 
ous anyway: "... if some other space than continuous functions is used for initial 
data, then the solution lies in C ...Therefore, for the fundamental theory, the 
space of initial data does not play a role which is too significant." [26, p. 33]. 
Anyway, Hale continues his comment by adding "However, in the applications, 
it is sometimes convenient to take initial data with fewer or more restrictions." 
In this sense, an alternative which has been quite studied is represented by 
the Hilbcrt product space C x L 2 (-r, 0; C d ), [3, 4, 16, 25, 35]. This second 
choice is often justified in the context of quadratic feedback control and linear 
filtering for retarded systems [15, 25, 38], for approximation reasons [28, 30], 
or when orthogonality is necessary [7]. In this manuscript the choice is that of 
continuous functions, reserving to present analogous arguments for the Hilbert 
space in forthcoming works of the authors, as already announced (and partially 
developed) in [11]. 

Once the proper state space is chosen, the long-time behavior of the evo- 
lution can be determined through the knowledge of the spectrum of infinite 
dimensional maps such as the semigroup of solution operators and its generator 
in the autonomous case, the monodromy operator for periodic problems, the 
evolution family in the nonautonomous case. The reduction of such operators 
to finite dimension allows to consider standard eigenvalue problems which can 
be easily solved, hopefully providing accurate estimates for the stability indica- 
tors (e.g. the rightmost root, the dominant multiplier or the largest exponent). 



x{t) = T{t,x t ), 



(1) 
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Usually, the construction of the finite dimensional approximation represents the 
"easy" step, although rather technical difficulties may arise when systems with 
possibly multiple discrete and distributed delays are the case. Instead, the the- 
oretical (i.e. not experimental) analysis of convergence is far away to be as 
simple and direct as the implementation of the numerical scheme. Here both 
aspects are equally treated with rigor. In particular, as far as convergence is 
concerned, the general properties that a underlying numerical method has to 
possess in order to lead to a satisfactory approximation of the spectral elements 
(eigenvalues, mutliplicities, eigenspaces, etc.) are recalled. To this aim, the 
(different) theories on spectral approximation of linear operators developed in 
the monographs [13, 24] are used as (alternative) background. 

Beyond this general treatment, the focus is placed on the pseudospectral 
approach, namely collocation together with polynomial interpolation: the core 
consists in substituting the exact operation (e.g. differentiation) to be done 
on a given function over a selected distribution of nodes with the same oper- 
ation as applied to the interpolating polynomial. The method, well-known for 
numerical integration of ordinary and partial differential equations as well as 
for the associated eigenvalue problems (see [36] for a guide), in the context of 
DDEs was first presented in [8] for the approximation of the spectrum of the 
infinitesimal generator for autonomous problems and in different guises also in 
[5, 9, 10, 12, 23, 37]. It benefits from the infinite regularity of the underlying 
eigenfunctions of the system at hand [26] and, by choosing to operate on Cheby- 
shev nodes, it is shown to be the unique one able to exploit all this regularity in 
problems with delay, too. As a result, the spectral elements are approximated 
at a convergence rate of infinite order, as it is rigorously demonstrated in the 
manuscript. To the best of the authors' knowledge, this represents the first (and 
only) complete proof of convergence in the field of DDEs. 

The paper is structured as follows. The prototypical model for (1) is in- 
troduced in Section 2 together with its state space representation through the 
associated evolution family. After some preliminaries discussed in Section 3, 
Section 4 deals with the numerical discretization. The convergence analysis is 
tackled in Section 5, precisely for the discretization scheme in Section 5.1 and 
for the spectral elements in Sections 5.2 and 5.3. In Section 6 applications are 
briefly discussed and commented. Appendix A collects Definitions and Lemmas 
useful and necessary to prove the main results of Section 5. 

2 Model, state space and evolution 

As a prototypical model for (1), we consider the scalar linear DDE with nonauto- 
nomous coefficients 

o 

x'{t) = a{t)x{t) + b(t)x(t - t) + [ c(t, 6)x{t + 6)d6 (2) 
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where r > is the maximum delay and [0, +00) 3 t i-> a(t),b(t) G C and 
[0, +00) x [— t, 0] 9 (i, 0) i-> c(i, 0) G C are continuous functions. All the 
arguments developed in the sequel apply as well to more general cases with 
matrix coefficients and multiple discrete or distributed delays, the extension 
concerning only technicalities useless to the treatment proposed in the paper. 

In order to focus on a well-posed Initial Value Problem (IVP) for (2), it is 
necessary to specify a suitable set of initial data. The classic literature [2, 18, 
26, 40] resorts to the Banach space of continuous functions (C, || • ||c) with 

C := C(-t, 0;C) 

and 

\W\c ~ g max. Q] M6)\- (3) 
Then, given two reals r > s and ip G C, the IVP 



x'(t) = a(t)x(t) + b(t)x(t - t) + J c(t, 6)x(t + 6)d6, t e [s, r] 
{x{s + 6)=y{6), 0e[-r,0] 



(4) 



admits the existence of a unique continuous solution x G [s — r, r] , which con- 
tinuously depends on the initial data [26]. This allows to introduce the linear 
and bounded operator T(r, s) : C — > C given by 

T(r, s)x s = x r , (5) 

x r , the state of the system at time r > s. In particular, x s — ip. The two- 
parameters family {T(r, s)} r > s is a strongly continuous evolution family and it 
is eventually compact, i.e. T(r, s) is compact as soon as r > s + t, see [26, 18] 
and [25] for definition and properties. 

In this work we are interested in the approximation of the spectral ele- 
ments of T(r, s), i.e. the nonzero eigenvalues as well as their multiplicities 
and eigenspaces, fundamental for stability purposes. 



3 Preliminaries and notation 

Depending on the role of the various mathematical objects, in general we use 
normal case for operators and functions (infinite dimension), bold case for ma- 
trices and vectors (finite dimension). 

After the time translation t^>-s + t, set 

C~ :=C = C(-t,0;C) 
C+ :=C(0,r 5 ;C) 
C* :=C(-T,r s ;C) 
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with r s := r — s (the same convention holds for function spaces other than C) 
and rewrite the IVP (4) as 

y'{t) = {G s y){t), te [0,r s ] 

y(0) = <p(0), 6e[-T,0], (b) 

for y(t) := x(s + t), x solution of (4), where the operator G s : — > C + is 
defined as 

o 

(G s y)(t) := a s (t)y(t) + b s (t)y(t - r) + J c s (t, 9)y(t + 6)dB, t £ [0, r s ], (7) 

— r 

where 

a s (t) := a(s + t), 
b s {t) :=b(s + t) 

and, for all 6 £ [-r,0], 

c s (t,6) :=c(s + t,6). 
The solution y £ in [— r, r s ] is intended as divided into 



y+(t):=y{t), t £ [0, 



:=y(t) = <p(t), t £ [-r,0]. 
For a given positive integer M, consider the grid of distinct nodes 

: = {~ r = : e M.,M <■■ < 0m,o '■= °i 

in [— t, 0] and set := C M+1 as the discrete counterpart of C~ , i.e. a function 
/~ £ C~ is discretized by the vector 

S M = ^-Af/ = (/ (^M,o)' • • • > / (^M,m)) T *= 

where : C~ — > is the restriction operator associating to a function its 
grid values at the nodes Sl^. Let : — > 11^ C C~ be the prolongation 
operator defined as 

M 

(V M v-)(t) :=^2iJt d (t)vJ, t £ [-t,0], 

for 

«- = (V,...,w M ) T £ C M 

and 
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the Lagrange basis polynomials relevant to the nodes . Then = V^jfM = 
^m^m/ - e is the polynomial of degree at most M interpolating /~ at 
the nodes O^, in fact 

7^ : C^ f — >• being the identity in , while 

^hi^-M = ^M' (8) 

£^ : C~ — > 11^ C C~ being the Lagrange interpolation operator on fl^. 
Similarly, for a given positive integer TV, let 

n+ :={0<8+ ,<■■■< 9+ N <r s } (9) 

be a grid of distinct nodes in (0, r s ), together with the auxiliary node 9~^ := 
(= 6^ ) and set := C N+1 as the discrete counterpart of C + , i.e. a function 
/+ e C + is discretized by the vector 

f + N = ^,o/ + = (/ + «o), ■ • ■ , / + «iv)) T e C+ 

where : C + — > is the restriction operator associating to a function 
its grid values at the nodes {8^ a } U Sl^. Let 7-*^ : — > C C + be the 
prolongation operator defined as 

M 

(V+ fi v+)(t) := ^ebjWvt, t e [0,r a ], 
i=o 

for 

»+ = («+..., 4) T eC+ 

and 

the Lagrange basis polynomials relevant to the nodes {9^ } U 0,^. Then 

/iv = ^Nfiftr = ^iv,o/ + ^/v,o G n jv is the polynomial of degree at most TV 
interpolating /+ at the nodes {9^ } U f2^, in fact 

1^ : — > Cff being the identity in , while 

V + 7? + — f + 

: C + — > 11^ C C + being the Lagrange interpolation operator on {9^ } U 
Six- In the sequel it will be necessary to refer rather to the Lagrange inter- 
polation operator on (i.e. without 9^ ). This latter will be denoted by 

: C + — > n^_ 1 C C + and given as 

£ + n = (io) 
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with obvious meaning of the restriction and prolongation operators relevant to 
fi+. 

Finally, for a given state y t € C, t G [0, r a ], and for the same integer M 
previously adopted for the grid fl M , we consider its discrete counterpart 

vt.u = (yt(6 M ,o)>--->yt( M,M^ T e Cm 

with Cm ■= C M+1 the discrete counterpart of the state space C. Correspond- 
ingly, 

M 
3=0 

is the relevant interpolating polynomial. It is not difficult to see that 

Vt,M = C M y t . 

Remark 1 Observe that C and C~ are the same space, as well as Cm and C M , 
all isomorphic to C M+1 . However, we reserve to distinguish the notation for 
the relevant meaning and role. Conversely, C~^ will be in general different from 
C M , since N ^ M can be chosen, their role becoming clear after the analysis of 
convergence in Section 5. 

4 Discretization of the semigroup 

We aim at finding a finite dimensional approximation T m,n{t, s) of the evo- 
lution family T(r, s) in (5). We basically use collocation, together with classic 
polynomial interpolation as introduced in Section 3. Briefly, in a discrete fashion 
we first transform the initial state into the solution on [0,r s ] and, second, we 
restrict this latter when r s > r, respectively prolong when r s < r. 

According to the notation set in Section 3 (but neglecting the reference to 
r s for simplicity) , we first construct matrices UJj M : C M — > C^ and N : 
C^ ->■ C^ such that 

Un,nPn = Ux,m<Pm (11) 
where pn,m G C ± is divided into 

fp+0) :=p N ,M(t), t G [0,r s ] - 12 - 

IPmW —Pjv.mW = <PM(i), t G [-r,0]. 

with determined by collocation of (6) on Q,^ with initial function ip M : 

(Pn)'KJ = ( G sPn,m){0+.), i = 1, . . . , N, 
p+(0) = ^(0). [6) 
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It is not difficult to check that the above matrices have entries, respectively, 
1 







9+ 

7 N,i 



{J: 



-*.KJt} f jK,i-'r) 





and 



[^JV,m]*.J : — 



1 





where 



'{ 

N+ = N+(r s ,T) := max {0+ . - r < 0} 

7 = 1,...,JV ' J 


















1,.. 


.,iV 


1 




n 

u, . . 




7V+ 


+ 1,.. 


n 

u, . . 


N 

















1,.. 


.,M 


1,.. 


.,JV+ 


0,.. 


.,M 


N+ + 1,. 


0,.. 


.,M, 



and is the Kronecker's delta. 

Second, and independently of the model coefficients a, b and c, we construct 
matrices V M N : — > Cm and V M M : C M — > Cm such that 



2/r,M — ^M.JvPjV + ^M,M(Cjlf 



(14) 



by restriction of pn,m to [r s — r, r s ] when r s > r, respectively prolongation by 
<Pm when r s < r. In particular, it is sufficient to define the above matrices with 
entries, respectively, 







i = 0,...,M~ 
j = 0,...,N 
i = M~ + 1, . 
j = 0,...,JV 
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and 



[V 



M,MlV ■- 




,M 



where 



M = M (r s ,r) := max {r s + 6 M > 0}, 

j—0,...,M ,J 



with the convention that V M N is full and V M M is empty when M~ = M, i.e. 
for r s > t. 

Eventually, by setting y Q M — ip M , it follows from (11) and (14) that 

Vt s m = T M ,N(r, s)y 0<M (15) 
is the sought discrete approximation of (5) with Tm n(t,s) : Cm — > Cm given 

by 

T M ,N(r,s) = V MN (U% tN ) U NM + V MM . 

Standard approximation arguments ensure that N is invertible for suffi- 
ciently large TV. 

Aim of this research is to show how and under which conditions the (com- 
putable) spectrum of Tm,n{i~, s) approximates that of T(r, s). 

Remark 2 Let us observe that suitable quadrature rules have to be applied 
whenever the integrals in the above matrices IT^ N and M cannot be com- 
puted exactly. Of course, there is no choice but that of uniform distribution in 
order to exploit the nodal values already at disposal due to the varying integration 
windows [9^ i — t, J, i = 1, . . . , N. However, we feel like to advice the use 
of gaussian-type nodes (e.g. Gauss-Legendre or Chebyshev): at the price of ex- 
tra technicalities in forming the matrix coefficients, this choice does not corrupt 
the spectral convergence that can be performed for the overall approximation as 
demonstrated in Section 5 (see [8] for further details). 



5 Convergence analysis 

The evolution family T(r, s) in (5) is an infinite dimensional map T(r, s) : C — > C, 
contrary to its matrix discretization Tm,n{t, s) : Cm — > Cm in (15). For compar- 
ison, it is therefore necessary to introduce an intermediate infinite dimensional, 
possibly finite rank, map Tjv(r, s) : C — > C. Set then 

T N (r,s)<p = (q N )r s (16) 

where € C ± is divided into 

/«£(*):=«*(*), *e[0,rj , , 

\q N (t) := q N (t) = ip{t), te[-r,0], k ' 
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with determined by collocation of (6) on Sl^ with initial function <p: 
Hti)'K.i) = (GsQN)K, i ), i = l,...,N, 

Note that, in general, the collocation polynomials gjy above and pn,m in (12) 
and (13) are different since relevant to different initial functions, ip and tp^, 
respectively. They coincide only when ip e 11^ . 

This Section is devoted first to provide in Section 5.1 error bounds for the 
remainder T(r, s) — Tjv(r, s) in a suitable state space. Such errors will be mea- 
sured in a pointwise sense in general (i.e. as applied to a given function in the 
chosen space), reserving to comment on the convergence in norm (i.e. as applied 
to all functions in the chosen space) . 

Second, in Section 5.2 a link for the spectral elements of TM,N{r,s) and 
IV (r, s) is studied, based on the relation 

T M , N {r, s) = TZ M T N (r, s)V M , (19) 

not difficult to be verified. 

Eventually, according to the theory developed in [13], it will be proved in 
Section 5.3 that pointwise convergence of Tjv(r, s) to T(r, s) in a suitable Banach 
space is a mandatory requirement for the approximation of the spectral elements 
of T(r, s) by those of Tjv(r, s), a finite number of which eventually coincide with 
those of Tm,n{t 1 s) as it will be demonstrated by virtue of (19). Instead, if the 
norm convergence is available (a much stringent requirement, rather difficult to 
happen), both the theories in [13] and [24] can be applied. These comment hold 
for any numerical method that can be potentially used to form the matrix ap- 
proximation TM,N(r 7 s) (e.g. Runge-Kutta or Linear Multistep based schemes, 
see e.g. [6, 20]). 

First to proceed, let us introduce some notation. For C and Y normed linear 
spaces, let B(C, Y) be the set of linear and bounded operators from C to Y. Far 
any A e B(C, Y) let 

II 4 II IWIk 

\\A\\c->Y ■= SUP m~7m ■ 

Vec \m\c 

If C = Y, then we will denote simply \\A\\c for the operator norm oi A G B(C) := 
B(C,C). 

For the convergence analysis, it will be often necessary to ask for more 
regularity than what so far demanded. To this aim, absolute continuity (see 
Definition 24 in Appendix A) will be used and the space 

C A := W 1 ' 1 = {tfj E L : $ £ L} 

will be considered where 

L := L\-t,0;C). 

Spaces C\ and C\ as well as L~ , L + and are similarly defined whether 
required. We will also resort to Lipschitz continuity, denoting the relevant spaces 
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with Lip^ and Lip^, Lip^- and Lip^- whether required (K denotes the Lipschitz 
constant). 

For C we use || • ||c as defined in (3), while for Ca we use 

H\\c a :=U\\l + W\\l (20) 

with 

o 

Ml = J \m\M- 

— T 

Similar norms for C~, C + and C^ 1 , C~^, C\ and C A and L~ , L + and L^ 1 are used 
whether required. With the above choices, all the spaces are of Banach type. 

A series of technical and preparatory Lemmas are stated and proved in 
Appendix A. Here we fully address only the main results. The constants C 
appearing in the various statements (included those in Appendix A), although 
different, always depend on r, s and r as well as on the coefficients a, b and c. 



5.1 Convergence of Tjv(r, s) to T(r, s) 

In this Section we study the convergence of the collocation under suitable hy- 
potheses on the coefficients a, b and c and on the initial function. 

Assumption 3 Assume the nodes in f2 J- to be the zeros of the N -degree Cheby- 
shev polynomial of the first kind in (0, r s ) [14], i.e. 

« .-^i-~fe^Yi.«-i s. 



'NA 



2N 



Assumption 4 Assume a s G Lip^, b s G Lip^ and c s (-,0) G Lip£ ,m for all 
9 G [— r,0] with K c G L. Assume, moreover, \\c s (- , 9)\\ c + G L and denote 

o 

\\Cs\\c+,L ■= J \\c s (;9)\\ c+ d9. 

— T 

Theorem 5 Let G s and be defined by (7) and (10), respectively, and let 
y G C ± be the solution of (6) with <p G C. Then, under Assumptions 3 and 4 
and for sufficiently large N, there exists a unique collocation solution qjy G 
given by (17) and (18). Moreover, 

\\y - q N \\e± < C\\p N \\ C+ (21) 

and 

\\y'-q' N \\c+ <C'\\ PN \\ C+ (22) 
hold where C and C are constants independent of N and ip and 

Pn ■= (I - C+)y' G C+. (23) 
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Proof. Given ip G C, let y G C ± be the solution of (6) on [— r, r s ]. It is not 
difficult to see, by integration, that y satisfies the functional equation in C ± 



y = u v + VG s y 
as soon as we consider u G as the function 



(24) 



MO), te [0,rj 



V : C + — > C ± as the integral operator 



(^)W ~ 4 



y y(cr)der, t G [0, r s ] 



(25) 



10, 



t G [-T, 0] 



and G s as given by (7). 

As for the collocation polynomial determined by (17) and (18), observe that 



q N (t)=q N (0) + 



J {<&)'(r)d° 



N 

¥>(<)) + jY,m+{o){q + N )'{0 + N .)do 



o i=i 



AT 

= ¥>(0) + / £mt(<7) \{G s q N )(6+j] da, t G [0,r s ], 
o i=i 

holds where to+'s are the Lagrange basis polynomials relevant to the nodes in 
fix (different from the ^ .'s). Then q N satisfies the functional equation in C ± 

qN = u v + VC^G s q N . (26) 

Consequently, the error function given by 

e N := y - q N (27) 

(which is clearly zero in [— r, 0]) satisfies the functional equation in C ± 

e N = VC+G s e N + Vp N , (28) 

as it can be seen by subtracting (26) from (24) and by adding and subtracting 
VC\jG s y in the result. 

Now, it is not difficult to see that (28) has a unique solution given by 



e^v = Ve 



N 



(29) 
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where e^v G C + is the unique solution of the functional equation in C + 

e N = C+G s Ve N + p N (30) 
thanks to Lemma 31. Moreover, 



e' N (t) :-- 



e N (t), t e [0,r s ] 

o, te[-T,o]. 



The thesis is now straightforward. ■ 

Observe that Theorem 5 shows that, under Assumptions 3 and 4 and for 
sufficiently large N, the continuity of the initial function ip is enough to guaran- 
tee the existence and uniqueness of the collocation solution as well as the error 
bound (21). However, from (23) it is clear that continuity is not sufficient to 
ensure convergence, i.e. 

lim \\y-q N \\ c ± =0, (31) 

N— s-oo 

for all p € C. In fact, G s y involves also the initial function <p due to the presence 
of the delay, hence it is only continuous and Faber's Theorem [14, 21] prevents 
convergence in all C ± . 

We now elaborate more, by resorting to the interpolation result in [33] , which 
holds only under Assumption 3. 

Proposition 6 Let pM be given by (23) under Assumption 3 and 4- If p> € Ca, 
then 

lim ||pjv||c+ = 0. 

N— >oc 

Proof. If ip G Ca, then G s y £ C\ follows easily from (7). The thesis is then 
given by the result in [33], valid under Assumption 3. ■ 

Proposition 6 shows that absolute continuity is a minimal assumption for 
the initial function ensuring convergence as meant in (31). 

Now we study the error for the approximated evolution family (16) w.r.t. 
the exact one (5). 

Corollary 7 Let T(r,s) be given by (5). Then, under Assumptions 3 and 4 
and for sufficiently large N, Tjy(r,s) in (16) is uniquely defined and, for any 
given <p G C, 

\\(T(r,s)-T N (r,s))<p\\ e < C\\p N \\ c+ 

holds where px is given by (23) and C is a constant independent of N and p. 
If in addition, p G Ca, then 

lim \\(T(r,s)-T N (r,s))cp\\c=0. 

N— >oo 

Proof. The thesis is straightforward by observing that 

\\(T(r,s) -T N (r,s))<p\\ c < \\y - qw\\c± 
holds for all r > s and by applying Theorem 5 and Proposition 6. ■ 
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Corrollary 7 ensures convergence for the approximated (still infinite dimen- 
sional) evolution family when applied to any absolutely continuous initial func- 
tion in the topology induced by || • ||c- In order to apply the theory developed 
in [13], the basic requirement to be satisfied is that of pointwise convergence in 
a Banach space. Precisely (and according to Definition 12 later on), given a 
Banach space (B, || • ||s), a sequence of operators {^4at}^ =0 m 1S pointwise 

convergent to A 6 B{B) if 



for all / 6 B. The property is denoted by An — > A. Clearly, Theorem 5 and the 
successive comment shows that T N (r,s) A T(r,s) cannot hold in the natural 
state space (C, || • \\c). We now prove that it holds instead if we choose as the 
state space (C A , \\ ■ \\c A ) with II • \\e A S iven b Y ( 20 )- 

Theorem 8 Let T(r,s) and TV(r, s) be given by (5) and (16), respectively, 
under Assumptions 3 and 4 and for sufficiently large N. Then Tjv(r, s) A 
T(r,s) in (Ca, \\ ■ \\c A )- Moreover, 



Proof. We use the same notation as in the proof of Theorem 5. Since 



pointwise convergence follows by virtue of Proposition 6. The last assertion 
follows from the Banach-Stcinhaus Theorem [13, Theorem 3.1]. ■ 

Let us underline again that the sequence of operators {Tpj(r, s)}^ =0 i s ma dc 
of infinite dimensional maps which have finite rank only when r s > t. In fact, 
when r s < r, the image Tjv(r, s)<p contains a piece of the initial function <p, 
precisely 



which, in general, is not a polynomial. However, for r s < r, T(r, s) is neither 
compact [26] . Therefore, for r s > r we have constructed a sequence of finite rank 
approximations to the exact evolution family, which can be proved to remain 
compact in B(Ca) as well by standard arguments. 

We conclude the Section with an important remark on the convergence in 
norm. 

Remark 9 Absolute continuity is sufficient to provide pointwise convergence, 
but for norm convergence more regularity is needed. In fact, Proposition 6 states 
that the interpolation remainder \\pn\\c+ vanishes, while for norm convergence 



Jim ||(A-^)/|| B = 



sup ||Tjv(r,s)|| Cj4 < oo. 



zn\\c a < r s (||ejv||c± + ll e w||c+) 




Pn\\c+ < g(N)\\y\\ c ± 



with g(N) vanishing independently of y would be necessary. 
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It is not difficult to see from the proof of Theorem 8 and from standard 
interpolation results as applied to that, if we further restrict to the state 
space C 1 := C 1 (— r, 0; C) with its natural norm 

U\y := iMic + ihft'Hc 

then we obtain convergence in norm, i.e. 

lim \\T(r,s)-T N (r,s)\\ c i = 0. 

N— >oo 

Further comments on this result will be made at the end of Section 5.3. 

5.2 Relation between T N (r,s) and T MjN (r,s) 

With an eye kept on the spectral elements, we analyze now the relation between 
the infinite dimensional approximation Tpj(r,s) defined in (16) and the finite 
dimensional one T m,n(t,s) introduced in Section 4 through (15). To this aim, 
N is considered fixed throughout the whole Section. The following result is 
fundamental to the scope. 

Proposition 10 Let the matrix Am '■ Cm Cm be given and define the oper- 
ator Am ■= V m Am71 m E B{C) through the restriction and prolongation maps 
IZ M and V M , respectively, introduced in Section 3. Then Am and Am have the 
same nonzero eigenvalues with the same geometric and partial multiplicities. 

Proof. Let \i e C \ {0}. We prove that 

V M ker(/iJ M - A M ) = kcr(/i/ - A M ). 

Let v £ Cm, b^O, such that 

Amv = /j,v. 

Let ip = V M v. Then 

A M ip = V M A M v = V M /J,v = flip. 
Vice versa, let ip e C, ip ^ 0, such that 

A M ip = f.ip- 

Since /i ^ 0, we have 

tp = V M v 

where 

v = -AmTI'm^ e c m- 

Hence 

V M A M v = nV M v 
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and then 

A M v = [iv. 

Now, we prove that 

dim ker (juTm — Am) = dimker(,u/ — Am) 

Let v\, . . . , v g be linearly independent elements of ker(/z/M — Am)- We have 
shown above that 

V'i = ^M^i G ker(/xJ - A M ), i = l,...,g. 

Since 

9 9 9 

=V M ^2 ai Vi =0=>^2aiVi = 0=> a l = 0, i=l,...,g, 

i—1 i—1 i—1 

the elements ipi, . . . , ip g are linearly independent. Viceversa, let tpi, . . . , tp g be 
linearly independent elements of ker(^7 — Am)- We have shown above that 

i>i = V M v u i = l,...,g, 

where Vi G kcr(filM — Am)- Since 

9 9 9 

i—1 2—1 z— 1 

the elements ui, . . . , v g are linearly independent. 

Finally, let us prove that there is a one-to-one correspondence between Jor- 
dan chains of Am and Am- Let vi, . . . ,v g be a Jordan chain of Am- Then 

— V M v i, • • • j "09 = 'Pm v 9 1S a Jordan chain for Am- In fact, ipi is an eigen- 
vector of Am and 

((il - A M )ipi+i = {nl - A M )V M Vi + \ 

= V M {fj,I M - A M )v i+ i 
= v t , i = 0,...,g-l. 

Vice versa, let Vi, . . . , ip g be a Jordan chain of Am- We have seen that 

for some «i G Cjv eigenvector of Am (since ipi is an eigenvector of Am)- Note 
that, for i = 0, . . . ,g — 1, if ^ = ? tf «j for some Vi G Cm, then ^> i+ i = P M v i+ i 
for some Uj+i G Cm- In fact, 

- A M )ipi+i =ipi= V M Vi 

and so 

i>i+i = V M ( -(i M K^i+i - «i) 
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We conclude that 

ipi = V M v t , i = l,...,g, 
for some Vi £ Cm, i-e. v\, . . . , v g is a Jordan chain for Am- In fact, 

T ? M V ^ =i>i = (v 1 - A M )^i+l 

= (fil - A M )V M v i+ i 

= V m {iiIm - A M )v i+ i, i = 0,...,g-l, 

and then 

Vi = (fiI M ~ A M )v i+ i, i = 0,...,g-l. 

■ 

Proposition 10 as applied to Tm,n{t, s) shows that this latter and the op- 
erator V M TM,N{r, s)lZ M have the same nonzero eigenvalues and multiplicities. 
Moreover, in the proof it can be seen how the eigenvectors of T m,n{t, s) and 
eigenfunctions of V M 1 'm,n{t, s)1Z m are related to each other via the restriction 
and prolongation operators. 

Now, observe that (19) and (8) give 

V M T M M{r,s)K M = C M T N (r,s)C M . 

In view of this and of the results of Section 5.1, we need only to investigate the 
relation between the spectrum of £ M T N (r, s)C M and that of T N (r, s). 

Theorem 11 Let Tjy(r,s) be defined by (16) under Assumptions 3 and 4 and 
with r s > t and let C M be defined by (8) with M > N. Then the spectral 
elements o/T/v(r, s) coincide with those of £ M TN(r 7 s)j£ M . 

Proof. On the one hand, we soon observe that for M > N and r s > r we have 

C M T N (r,s)=T N (r,s). 

In fact, for r s > r the range of T N (r, s) is 11^ and £^11^ = 11^ for M > N. 

On the other hand, let ip be an eigenfunction of T N (r,s) associated to the 
eigenvalue /x, i.e. 

T N (r 7 s)tp = flip. 
Then, for r s > t, it must be t[) £ 11^. Consequently, 

T N (r, s)£ M ip = T N (r, s)ip 

whenever M > N. Arguments similar to those used to prove Proposition 10 
complete the proof. ■ 

As an immediate consequence of Theorem 11, it is sufficient to take M = N 
in the construction of the matrix Tm,jv(t*, s) in Section 4, so to keep the compu- 
tational effort as low as possible. Moreover, it is interesting to observe that no 
special assumption on the distribution of the nodes VL M has to be made, con- 
trary to Assumption 3 for f2^. In fact, these nodes serve merely for representing 
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Tj\r(r, s) in finite dimension, and not for approximation reasons as those in fi^. 
However, Theorem 11 holds only for r s > r, which is the case of interest since, 
as already remarked, this condition guarantees the compactness of T(r, s), and 
hence the nice and known properties of its spectrum. Nevertheless, it is not 
difficult to see that pointwise convergence £^Tjv(r, s)C^j A T/v(r, s) holds in 
all (Ca, II • \\c A ) as M — > oo and N fixed as long as r s > r. 

5.3 Computation and convergence for the spectrum of 

T(r, s) 

From now on we assume r s > r and M > N. 

In Section 5.2 it is shown that, for a fixed N, the spectral elements of Tjv(r, s) 
coincide with those of the matrix Tm,n(x, s), through which they can be effec- 
tively and efficiently computed as the result of a standard algebraic eigenvalue 
problem in finite dimension. The last step consists in proving that these ele- 
ments converge to those of T(r, s) in the limit as N —> oo. To this aim we apply 
the theory developed in [13], that we briefly recall for what concerns its basic 
facts and as adapted to the notation of the present manuscript. 

Following [13], let fi be a nonzero and isolated eigenvalue of T(r,s) with 
finite algebraic multiplicity m, geometric multiplicity g and ascent I. Let A be 
a neighborhood of \x such that \i is the only eigenvalue of T(r, s) in it. Let T be 
a closed Jordan curve isolating fi and drawn in A. Observe that all this makes 
sense since, for r s > r, T(r, s) is compact and hence it has only point spectrum 
with nontrivial isolated eigenvalues of finite algebraic multiplicity plus possibly 
as accumulation point [13, Theorem 2.34], [26, Chapter 7]. To such a /i it is 
associated the spectral projection 

P:=^:£(T(r,s)-zI)- 1 dz 

and the relevant generalized eigenspacc A4 := PC a- Recall that m = dimM, 
g = dimker(T(r, s) — fil), while (T(r, s) — /x/)P is nilpotent and I is the maxi- 
mum integer such that ((T(r, s) — fiI)P) e = 0. Moreover, 1 < g, I < m < oo. If 
Q is another porjection, then for Af := QCa it is defined the gap 

B(M,N) := max{dist(M,AA),dist(A^,7W)} 

with 

<Kst(M,N) := sup dist(^,7V)^ sup dist(£,.M) =: dist(Af,M). 

\W\c A =i K\\c A =i 

The following definitions are necessary, also to understand the role of the 
analysis so far conducted. 

Definition 12 (pointwise convergence [13]) T^ir, s) A T(r,s) iff, for all 
f e Ca, Tfj(r, s)ip — > T(r, s)ip as N — > oo. 
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Definition 13 (stable convergence [13]) T/v(r, s) — zl A T(r,s) — zl for 

z e A \ Ox} iff 

(i) T N (r,s) A T(r,s); 

(ii) 3K > 0, 3N s.t. for N >N, (T N {r,s) - zl)- 1 G B(C A ) and \\(T N (r,s) - 

Definition 14 (strong stable convergence [13]) Tjv(r, s) — zl -—> T(r, s) — 
zl /or z e A \ {^} i/f 

(i) T N (r, 8)-zI±> T(r, s) - zl for z G A \ {/x}; 

(mJ dim PnCa = Tn for N large enough. 

IfTjy(r 7 s) — zl T(r,s) — zl in A, then Tjy(r,s) is said a strongly stable 
approximation ofT(r,s) in A. 

Definitions 13 and 14 above hold similarly for all z G T if both conditions 
(i) hold for z G T. 

A first important result follows. 

Proposition 15 (Proposition 5.6 in [13]) 7/TV(r, s) is a strongly stable ap- 
proximation ofT(r, s) in A, then, Tpj(r, s) has in A, for N large enough, exactly 
m eigenvalues, counting their multiplicities. 

Let us call [ijn, j = 1, . . . , m, such eigenvalues, let /in represent any among 
the distinct ones and let ipN, H^jvIIca = 1; be an associated eigenfunction. For 
this [j,n, let P/v and A4n be, respectively, the relevant spectral projection and 
generalized eigenspace as previously introduced. Since T(r, s) is not self-adjoint 
in general, /x is best approximated by the arithmetic mean 

m 

flN '■= ^ fJ-jN- 
= 1 

We recall now the two fundamental results, which involve the quantity 

e N := \\(T(r,s)-T N (r,s))P\\c A . (32) 

Observe that En is the remainder of the approximation as restricted to the 
generalized eigenspace, rather than on the whole space. 

Theorem 16 (Theorem 6.6 in [13]) IfT]y(r,s) is a strongly stable approxi- 
mation of T(r, s) in A, then, for N large enough, the quantities 

- P)1>n\\c a , for ipN G M N , Hn\\c a = 1, 
W - Pn)1>\\c a , foryjeM, 
Q(M,M N ), 

fj, - /X;v, 

it m n jN J 
are at least of order en- 
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Theorem 17 (Theorem 6.7 in [13]) //Tjv(r, s) is a strongly stable approxi- 
mation ofT(r,s) in T, then, for N large enough, 

max \fi - fi jN \ = 0(s]{ e ), 

j=l,...,m 

min \n-fi jN \ = 0(e 9 J m ), 

] = !,.. .,m 

dist(Viv,ker(T(r,s) - fxl)) = 0(s]( e ). 

In order for the previous results to hold true, we need now to verify the (only) 
hypothesis of strongly stable convergence according to Definition 14. This, in 
turn, requires (ii) in Definition 14, (ii) in Definition 13 and pointwise conver- 
gence. 

The latter is provided by Theorem 8, justifying all the analysis performed 
in Section 5.1. 

As for (ii) in Definition 13, it is enough to observe that z e A \ {^} implies 
z 6 p(Tjy(r, s)), the resolvent set of Tw(r, s), hence (ii) holds by the definition 
of this latter. 

Eventually, as for (ii) in Definition 14, we first recall the following definitions 
and results about the spectral projections Pn and P. Below we use B := {tp e 

c A ■■ \m\c A < i}. 

Definition 18 (collectively compact convergence [13]) Pn ^ P iff 

(i) P N A P; 

oo 

(ii) the set U (P — Pn)B is relatively compact in Ca- 
n=i 

Definition 19 (compact convergence [13]) P N A P iff 

(i) P N A P ; 

(ii) for any sequence {£jv}jv = i in B, the sequence {(P — Pn)£n}n=i * s rela- 
tively compact in Ca- 

Theorem 20 (Theorem 3.9 in [32]) The projections P and Pn are com- 
pact. 

Proposition 21 ([1], Proposition 3.13 in [13]) For projections P and Pn 
such that P is compact, the following are equivalent: 

(i) Pn A P and dim PnCa = dim PC a < oo for N large enough; 

(n) Pn ^ P- 

Theorem 20 ensures that Pn — P is compact, hence Definitions 18 and 19 
are equivalent, [13, p. 125]. Moreover, (ii) in Definition 19 is trivially satisfied. 
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Then, by virtue of Proposition 21, (ii) in Definition 14 is guaranteed as soon as 
Pjv A P. As for this latter, for any given ip £ Ca, 



and hence Pm — > -P follows from T/v(r, s) — > T(r, s), again guaranteed by The- 
orem 8. 

Now we are able to comment on the convergence rate. Thanks to Theorems 
16 and 17, the error between the spectral elements of Tjy(r,s) and T(r, s) de- 
creases as fast as (32). Thanks to Theorem 5, this latter is governed by ||/9jv||c+ 
with pm given by (23) for y solution of (6) with a generalized eigenfunction 
associated to /i as initial function. As it is well-known, the eigenfunctions of 
T(r, s) are analytic [26]. This, together with Jackson's type Theorems, ensures 
spectral accuracy, i.e. a convergence of infinite order: 



where K = K(\/i\) is a constant independent of N. 

We resume all the spectral convergence analysis in the following, where So(0 
denotes formally the spectral elements of an operator (i.e. any of the quantities 
mentioned in Theorems 16 or 17). 

Theorem 22 Let T(r, s) be given by (5) with r s > t. Under Assumption 3 and 
for sufficiently large N, T/v(r, s) in (16) is uniquely defined. For M > N , set 
T m ,n{t, s) = K M T N (r, s)V M - Then 

(i) for a fixed N, £ (?M,Jv(r, s)) = S (T w (r,s)); 

(ii) as N — V oo, So(P/v( r j s )) S (T(r, s)) with spectral accuracy. 

Let us underline that the equality stated in (i) above is true for the eigenval- 
ues and relevant multiplicities, while it hides the restriction and prolongation 
operators as far as the eigenvectors and eigenfunctions are concerned. 

Remark 23 By recalling Remark 9, the result about convergence in norm there 
stated may be used to apply the theory developed in [24] in order to obtain the 
similar results as given in this Section for the approximation of the spectral 
elements. Such a convergence analysis is similar to that developed in [10] for 
partial retarded functional differential equations, and it is not based on the theory 
developed in [13]. However, it is worthy to underline that we chose to follow 
[13] since it requires the least possible restriction of the natural state space, i.e. 
Ca rather than C . 
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6 Applications 



It is well-known (e.g. [26, 18]) that the asymptotic stability of stationary solu- 
tions can be characterized through the knowledge of the spectrum of T(r, s). In 
particular, the following two situations are certainly worthy to be mentioned in 
view of their importance in a wide class of applications. 

• When (2) is the result of the linearization of a nonlinear system of DDEs 
around an equilibrium solution, then the coefficients are autonomous, i.e. 
for all t e [0, r], a(t) = a, b(t) = b and c(t, 9) = c{6) for all 9 G [-r, 0] and, 
consequently, the evolution family {T(r, s)} r > s reduces to the standard 
Co-semigroup of solution operators {T(r)} r > . Moreover, for r > t, such 
an equilibria is asymptotically stable iff < 1 for all the eigenvalues \i 
of T(r), the so-called multipliers. 

• When (2) is the result of the linearization of a nonlinear system of DDEs 
around a limit cycle with period oj, then the coefficients are periodic with 
the same period, i.e. for all t e [0,r], a(t + to) = a(t), b(t + u) = b(i) 
and c{t + u,9) = c(t, 9) for all 9 e [— r, 0] and, consequently, the spectral 
properties of the evolution family {T(r, s)} r > s can be studied through the 
monodromy operator U{uo) := T(uj,0). Moreover, for any lo > t, such a 
periodic orbit is asymptotically stable iff < 1 for all the eigenvalues n 
of U(cu), the so-called Floquet multipliers. 

It is then clear that the results of the present work can be straightforwardly 
used in both cases without any additional effort. Let us observe, moreover, 
that for periodic problems where the period of the coefficients lo is less then 
the maximum delay r the monodromy operator is not compact, even though a 
sufficiently large power k of it is so, precisely U k (u) = U(kuj) with some integer 
k such that kuj > t. Then, for stability purposes, spectral approximations of 
the latter are sufficient. 

A Appendix 

Definition 24 Let (C,d) be a metric space. A function E D I 3 t i-> f(t) E C 
is absolutely continuous if for any e > there exists a 5 = 5{e) > such that 
for any sequence {[a„, fln]}%Lo of pair wise disjoint subintervals of I satisfying 




it follows 



oo 



Ew./W)< e - 



Lemma 25 Let V be defined by (25). Then 



\\V\\ e+ ^c± =r, 
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Proof. The inequality 

\\V\\ c +^c± <r s 

is trivial. For the equality it is enough to use the constant function 1 G C + . 
Lemma 26 Let V be defined by (25). Then Vy G for all y G C+ . 

Proof. For y G C + define, according to (25), z G as 

t 

z(t) := (Vz)(t) = { 



J y(cr)d<T, t G [0,rj 
t G [-r, 0] 



o 

10, 



By following Definition 24, let {[ot ni fi n }}^Lo be any sequence of pairwise disjoint 
subintervals of [— r, r s ] satisfying 

oo 

^2\(3 n -a n \ <6 
for a given 5 > 0. For a fixed index n it holds 



z(/3 n ) - z(a n ) = < 



r if a n < fi n < 

J y(a)dcr if a„ < < ^„ 
o 

y y(cr)da if < a„ < /?„. 



Then it easily follows 

\z(/3 n ) - z(a n )\ < \/3 n - a n \\\y\\ c + 
independently of whether [a n ,/3 n ] falls into [— r, r s ]. Consequently, 

oo oo 

- a n\\\y\\c+ < s \\y\\c+ 

n=0 n=0 

and hence y satisfies Definition 24 by setting 

e 



(33) 



\y\\c+ 



for any e > 0. ■ 

Lemma 27 Let G s and V be defined by (7) and (25), respectively. Then, under 
Assumption 4, G s Vy G C\ for all y G C + . 
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Proof. For y E C + define, according to (7), z E C + as 
z{t) := (G s Vy)(t)=a s (t)(Vy)(t) 

o 

+b s (t)(Vy)(t-T) + J c s (t,e)(Vy)(t + 6)d6, ie[0,rj. 

— T 

Then, according to (25), we have 

o 

z(t) = a s (t)(Vy){t) + J c s {t,6){Vy){t + 6)d6, te[0,r a ], 
-t 

whenever r s < t while 



z(t) = a s {t){Vy){t) + < 



U 

b s (t)(Vy)(t - t) + J c s (t,9)(Vy)(t + 8)d6, t E [r, r s ] 

— r 



J c s (t,9)(Vy)(t + 9)d9, t€[0,T] 



whenever r s > t. By following Definition 24, let {[«„, /3n]}^ any sequence 
of pairwise disjoint subintervals of [0, r s ] satisfying 



\Pn ~ a n \ < S 

for a given S > 0. Clearly, 

oo 

J2\(Pn + 0)-(a n + 9)\<5 

holds for all 9 E [— r, 0]. It is not difficult to see that, for a fixed index n, the 
bound 

\z(p n ) - z(a n )\ < |a s (/3„)| \(Vy)((3 n ) - (Vy)(a n )\ 
+ \a s (/3 n ) - a s (a n )| \(Vy)(a n )\ 
+|6,(/8n)||(V r i/)( j 8n-T)-(Vi/)(a n -r)| 
+\b s {p n )-b s {a n )\\{Vy){a n -T)\ 
o 

\c s (p n ,9)\ \(Vy)(p n + 9)- (Vy)(a n + 9)\ d9 



+ 



J \c s (p n ,9) - c s (a n ,9)\ \(Vy)(a n + 9)\ d9 
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holds thanks to (25) both for r s < t and r s > r. Now, by applying Lemma 25, 
(33) in Lemma 26 and by using Assumption 4 we get 

oo 

J2\*(Pn) - z(a n )\ < Sk\\y\\ c+ 



n=0 



with 



k ■■= [(Klle+ + ||Mc+ + \\c s \\c+,L)+r s (K a + K b + \\K c \\ L )]. (34) 
Hence z satisfies Definition 24 by setting 

S 



k\\v\ 



c+ 



for any e > 0. 



Lemma 28 Let G s , V and k be defined by (7), (25) and (34), respectively. 
Then, under Assumption 4, G s Vy <E LipJ- for all y E C + with K < k\\y\\ c + ■ 

Proof. For y 6 C + , Lemma 27 ensures that G s Vy is absolutely continuos in 
[0,r s ], hence also a.e. differentiable. Then it is sufficient to show that G s Vy 
has a bounded first derivative whenever defined. To this aim observe that for 
z = G s Vy defined as in the proof of Lemma 27 



z'(t) = a' s (t)(Vy)(t) + a s (t)(Vy)\t) + c s (t, -t)(Vy)(0) 
o 

^(t, 9){Vy)(t + 9) + c s (t, 9){Vy)'{t + 9) 



+ 



U 

/ 



de 



holds for a. a. t g [0,r s ] whenever r s < t. In fact, under Assumption 4, all coef- 
ficients are Lipschitz continuous, hence a.e. differentiable with bounded deriva- 
tive. Similarly, 



z'(t)=a' s {t){Vy){t) + a s (t)(Vy)'(t) + 



( o 

/ 



b' s (t)(Vy)(t-T) + b s (t)(Vyy(t-r) 
c s (t,-t)(Vy)(0) 



+ < 



dc 6 
dt 



dc s 
dt 



(t,9)(Vy)(t + 9)+c s {t,9){Vy)'(t + 9) 



(t,9)(Vy)(t + 9) + c s {t,9)(Vy)'(t + 9) 



d9, in [t, r s ] 
d9, in [0, t] 



holds for a. a. t G [0,r s ] whenever r s > r and the derivatives of the coefficients 
are bounded again. Since 

(VyY(t) = y (t) 

for all t E [0,r s ] and 

(Vy)(0) = 0, 
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it follows 



z'(t)=a' s (t)(Vy)(t) + a s (t)y(t) 
o 

/ ^(t,e)(Vy)(t + 6)+c s (t,6)y(t + 0) 



for a.a. t <E [0, r s ] whenever r s < r while 



<W 



z'(t)=a' s (t)(Vy)(t) + a s (t)y(t) 
o 

-^(t,6)(Vy)(t + 6)+c s (t,e)y(t + e) 



U 

+/ 

— T 



d(9 



(t)(Vy)(t - r) + b s (t)y(t - r), in [r, r s ] 

in [0,t] 



for a.a. t <G [0, r s ] whenever r s > r. This shows that z' is wcll-dchned in all 
[0, r s ] except for the set of measure zero where the derivative of the coefficients 
is not defined and also for t = r when r s > r. In fact 

z'(t+)-z'(t-) = 6 s (t)2/(0) 

and, in general, y(0) ^ 0. Nevertheless, for all the other values we have 

\z'(t)\<k\\y\\ c + 

with k in (34). Since the same bound is valid also for r s < r, the proof is 
complete. ■ 

Lemma 29 Let G s , V and be defined by (7), (25) and (10), respectively. 
Then, under Assumption 4, 



\\(I-£+)G s V\\ c+ <C^ (So) 



where C is a constant independent of N and A^- is the Lebesgue constant relevant 
to the nodes in (9). If, in addition, Assumption 3 holds, then 

lim \\(I-£+)G s V\\ c+ =0. (36) 

N— >co 

Proof. By standard interpolation results such as Jackson's Theorem [14, The- 
orem 13.3.7], for any given y e C + , we have 



||(7 - £+ )G s Vy\\ c+ < (1 + k+)E N (G s Vy) 



<(l + k + N )[l + \\u(G s V y] ^-\ 



2N J 
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where En is the best uniform approximation error of G s Vy in IT^ and u>, 
the modulus of continuity, is bounded by virtue of [14, Theorem 1.5.1] and by 
applying Lemma 28. If, in addition, Assumption 3 holds, then (36) follows by 
observing that 

A+ =0(logJV) 
by virtue of Natanson's Theorem [34]. ■ 

Lemma 30 Let G s and V be defined by (7) and (25), respectively. Then, under 
Assumption 4, (7 - GgV)^ 1 e B(C+). 

Proof. We have to prove that for any given f E C + there exists a unique y G C + 
solution of 

(I - G s V)y = f. 

By contradiction suppose that, given /, there exist x ^ y such that 
(I - G s V)y = / = (/- G s V)x 

and set z := y — x. 

If r s < t, then z solves the Volterra Integral Equation (VIE) 

z(t) - J k s (t,a)z(a)da = 0, t € [0,r s ], (37) 



with 



u 

k s (t,a)=a s (t)+ J c s (t,9)d9, 



-t+a 

as it can be verified by a standard change of integration order for the distributed 
delay term. The VIE has a continuous kernel k s , hence [32, Theorem 3.12] 
ensures z = 0. 

If r s > t, then z solves the VIE 



t 



z(t)- [k s (t,a)z(a)da-\ J ^ + J c s (t,0)dO z(a)da, t E [r,r s ] 

J V -t+a ) 

° U t€[0,T]. 

(38) 

This latter does not have a continuous kernel due to the presence of the discrete 
delay term, hence [32, Theorem 3.12] cannot be applied. Nevertheless we can 
adopt a step-by-step procedure. In fact, being r s > r, there exist an integer 
a and a f3 6 [0,r) such that r s = ar + (3. We prove that z(t) = for all 
t € [(i — 1)t, ir], i = 1, . . . , a, and also for t e [or, r] . For i = 1 we have i € [0, r] 
and the VIE (38) reduces to (37) for which [32, Theorem 3.12] provides z(t) = 
for all t e [0, t}. For i = 2we have t S [r, 2r] and the last integral in (38) over 
[0, t — t] involves the solution over [0, t] , which has just been proven to be null, 
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showing that z(i) — for all t E [r, 2r] as well. Iteration of the process leads 
to z(t) = for all t E [0,r s ] since the same argument holds for the last interval 
[ar, r], too. 

Eventually, independently of whether r s < r or r s > r, we get z = and 
y = x, which contradicts the initial assumption. ■ 

Lemma 31 Let G s , V and be defined by (7), (25) and (10), respectively. 
Then, under Assumptions 3 and 4 and for sufficiently large N , {I—C^GgV)^ 1 E 
B{C + ). In particular, 

\\(I - C+GsV)- 1 ^ <2\\(I-G S V)- 1 \\ C+ . 

Proof. The thesis follows immediately by observing that 

I - C+G S V = (I - G S V) + (I - C+)G S V 

and by applying the Banach's perturbation Lemma [32, Theorem 10.1] whose 
hypotheses are satisfied thanks to Lemmas 29 and 30. ■ 
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